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Abstract 

We calculate the rates and lepton (£) invariant mass distributions for decays 
of the form ^ — > £ + £~j, which are important radiative corrections to the 
purely leptonic decays h — ► £ + £~. Our approach uses the loop diagrams 
which arise by including the two photon intermediate state and we retain the 
imaginary parts of the loops - a radiative extension of the 'unitarity bound' 
for the process. These results are compared which those obtained using a 
model in which the meson couples directly to the leptons. 
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I. INTRODUCTION 



Recently, a few electron events in the channel K\ — > e + e~, where is the long lived 
neutral K meson, have been observed The same experiment sees thousands of muon 
events, K\ — > One of the experimental acceptance conditions is that the invariant 

mass of the two leptons be within a few MeV of the K mass. This means events are lost if 
a photon of sufficient energy is also emitted, i.e. K\ — > £ + £~j, with £ = e or fi. 

Our purpose is to estimate the size of this radiative correction by calculating 

1 am— a) 2 dr 

— / ds—. 1 

T Jim? ds 

Here V is the rate for — > £ + £~ , ~f, T is the rate for K\ —>■ £ + £~ and s in the square of 
the lepton invariant mass, s = (p + p') 2 , where p and p' are the momenta of the lepton and 
antilepton. The K mass is denoted by M, the lepton mass by m and A is an experimental 
parameter. In the next Section, we examine a tree level model which treats the meson-lepton 
interaction as a point coupling. Section III contains results for a model in which the meson- 
lepton point interaction is replaced by a loop diagram with two photons in the intermediate 
state. We conclude with a comparison of the two approaches. Corrections to K\ — > 7r + 7r~7 
are presented in an Appendix. 



II. RADIATIVE CORRECTIONS IN THE TREE APPROXIMATION 

The differential decay width dT /ds has been calculated in a simple model @ where the 
meson-lepton coupling is taken to be a pseudoscalar interaction with an effective coupling 
constant g. This leads to a transition amplitude T of the form 

T = 9 e ("M / m ^ +y W o + « W m y W )) , (2) 

as illustrated in Fig. 1. Here, e is the photon polarization vector and k is the photon 
momentum. Unfortunately the result is given only in the limit M 3> m, a condition that is 
clearly not satisfied for K decay into muons. The extension to include terms of all orders 
in m 2 /M 2 is straightforward and we find that Eq. (0) gives a differential decay width for 
K° L -> £ + £~~i of the form 



1 dT a 1 1 



T ds ttM 4 ^\-Am 2 /M 2 M 2 - s 



(M 4 + s 2 - AM 2 m 2 ) In ( ) - 2M 2 sv 

1 — v 



(3) 



where v = J 1 — 4m 2 / s. This is a slight extension of Bergstrom's expression || by the terms 
proportional to m in the numerator and denominator. 

In an attempt to estimate the model dependence of our corrections, we will compare 
our results with "model independent" corrections given by keeping only the universal soft 
bremsstrahlung correction terms of the form 

T - 9 eu(p) l5 v(p') fef " ^) ■ (4) 
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In this case the contents of the square brackets in Eq. (|3|) are replaced by 

'1 + 



(2s 2 - Asm 2 ) In 



1-v 



2vs 2 . 



(5) 



For the expressions (|3|) and ([5]), the integral in Eq. ([[]) can be evaluated analytically to 
determine the fraction of lepton pairs missed. Using Eq. @ we get 



1 f(M—A.) 2 dT a 
as ■ 



Tq J Am 2 



where 



ds n a/1 - e- 



2A A 2 _2m 

17 IT 2 ' 17 



F(5,e) 



(6) 



and 



F(S, e) = [-2(2 - e 2 ) (in 6 + 2 ln(l + (3 2 _)) - 2 (l - 5 + |(1 - <5) 2 ) + 2(1 - VT^ 2 ) 2 



+e 2 (l + le 2 )] lna+ + (2- £ 2 



_ Li2 K + Li2 « 



an 



-2 Vl - £ 2 In 



with 



l-a%./P: 

Vl-5 y/l-5-e 
a + = : 1 : 



(7) 



e e 



(8) 
(9) 



Li 2 (x) is a Spence function or dilogarithm defined as 

r x Ht 

Li 2 (x) = - -ln(l-t). 
Jo t 

Numerical values of Li 2 (a;) can easily be obtained using Maple or Mathematica. 
For expression (|5|) we again have Eq. (^) where now 

F pole (5, e) = [-2(2 - e 2 ) (in 5 + 2 ln(l + /3 2 )) - 2 ((2 - e 2 ){\ - 5) + (1 - 5) 1 



+2(1 - v / T^7 2 ) 2 + £ 2 (2 - | £ 2 )j In a + + (2 - e 



-Li 2 



+ Li 2 



an 



-2V1 -e 2 ln 



le 2 )VT^5Vl-5-e 2 . 



(10) 
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III. RADIATIVE CORRECTIONS USING A ONE-LOOP MODEL 



All the experimental results announced so far M find the rate for K\ — > fi + fi to be 
near the theoretical lower limit given by multiplying the rate for K\ — > 77 by the rate for 
77 -> fJ- + »~ 



r = r(K^77)^ r 




(ii) 



where 



P = ^/l-4m 2 /M 2 . 

The rate for K\ — > e + e~ as given in U is larger than the unitarity bound given by Eq. flTTD 
and consistent with predictions from chiral perturbation theory j3J. Nevertheless, it seems 
reasonable, in attempting to extend the calculation of Eq. (P beyond the result obtained 
using Eq. (fj), to calculate the absorptive part of K\ — > £ + £~ / y diagrams shown in Fig. 2. 
In particular the box diagram takes us beyond simple bremsstrahlung off external legs. 
Like Eq. ([TT|), the bremsstrahlung pole terms of Fig. (2), which vary as l/oo, where u is the 
photon energy, contain a factor of the lepton mass. It has been suggested that the terms 
which vary as uj to a positive power might not include a lepton mass factor and could thus 
be anomalously large ||. 

To evaluate the diagrams of Fig. 2, we take the if -photon-photon effective Lagrangian 
to be 

Ay 7 F^Ffa, , (12) 

where F^ v is the photon field tensor, F^ u is its dual, is the K\ field and Ay 7 is a constant. 
This leads to the vertex function 

r^ t! /(/c, k ) = 2A 77 £[iuaf3 k k , (13) 

where k and k! are the photon momenta. In general, this expression can also include a form 
factor which depends on k 2 and k' 2 . The implications of including this additional factor are 
discussed below. 

In this case, the expression for dT/ds is very complicated and we will not attempt to 
write it out. The box graph involves integrals of the form 

d A q t(l v t , (f(l v , (f / 14 x 

q 2 [(q + p) 2 — m 2 ] [(q + p + k) 2 — m 2 ] (q + P) 2 

The triangle graphs involve similar integrals with one or two factors of q in the numerator and 
either the second or third factor in the denominator of Eq. ( JH| ) omitted. These integrals can 
be expanded in terms of the external momenta as outlined in the appendix of Passarino and 
Veltman ||. The momentum expansion and its scalar coefficients are given by a computer 
code called LOOP []7||| which is a slight modification of a code written by Veltman called 
FORM Factor. Within this code the integrals are evaluated in terms of Spence functions as 
defined above and these functions are then evaluated numerically. 
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Once the amplitudes are determined they are squared and summed over spin, including 
the photon polarization, in the usual way. To check for errors we replace the photon polar- 
ization by its momentum and look for gauge invar iance. This is our only real check but it is 
a very powerful one since gauge invariance requires a delicate cancellation among the three 
diagrams 0. 

The real part of the amplitudes diverges because the effective coupling Eq. (|T2|) has 
dimension 5. The absorptive part has several contributions: K — > 77 followed by 77 — > 
£ + £"7 as well as K — > •y£ + £~ followed by — > j£. This is illustrated by the cut diagrams of 
Fig. 3. In the first of these diagrams, the intermediate photons are on-shell (k 2 = k' 2 = 0), 
which is equivalent to our assumption that A 77 is constant. In the second, one of the 
photons is virtual, and the effective coupling has the general form A 77 » = A yy f(k 2 ), where 
f{k 2 ) is a form factor normalized to f(0) = 1 [Kj. Our numerical calculation of the complete 
absorptive part cannot separate these contributions, so we have effectively assumed A 77 . = 
A 77 throughout. This assumption is justified in the case of electrons, since the form factor 
correction to the width of the Dalitz decay K\ — + e + e~7 is only a few percent due to the 
preference for low e + e~ invariant mass. For the muon case, the form factor correction to the 
Dalitz decay width is 20-25% fflOf] , and there could be a discernable effect in the one loop 
contributions to dT/ds. 



IV. RESULTS AND CONCLUSIONS 

In Figs. 4, 5 and 6, we show the differential width obtained from the absorptive part of 
the one-loop calculation as a function of the lepton invariant mass, s, for K® L — > e + e~7, 
K° L — * /i + yU~7, and tt° — > e + e~ / y. For comparison we also plot Eqs. (Q) and (|^). The result 
of the loop calculation is almost the same as Bergstrom's differential width, as modified by 
us to include the lepton mass, and both differ substantially from the "model independent" 
width where only the 1/u terms are kept. 

In Tables I, II, and III we give the integrated width, Eq. (1), for several values of A. Again 
the result from the loop calculation is very similar to that given by Eq. (^) and quite different 
from that given by Eq. fllCf). For electrons the correction is large, but not anomalously so, 
and there is no indication that the diagrams of Fig. 2 are not proportional to the lepton 
mass. 

The moral would seem to be that, except for very small invariant masses, the more 
complicated calculation of the loop model of Sec. Ill is unnecessary and Eqs. (^) and (||) 
are sufficient. We have made no attempt to calculate the radiative corrections within the 
acceptance bin, (M — A) 2 < s < M 2 . For the model of Sec. II Bergstrom has given a 
complete expression for the correction from virtual photons. This, together with Eq. @, is 
all that is needed. To calculate the virtual corrections for the model of Sec. Ill is beyond 
the scope of this work. 
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APPENDIX: 



In a model for the K° — > 7r + 7r - vertex similar to Eq. (Q), where the Ktttt vertex is taken 
as a constant, the differential rate for K° — > 7r + 7r~7 has only the 1/u terms and is therefore 
given by Eq. (J|) multiplied by M 2 / s to remove the spinor factor | TT| . The integrated rate is 
given by Eq. (0) with 



-Ppoie(<5, e) 



+2e 



2(2 - e 2 ) (ln5 + 2 ln(l + 0*_)) - 4(1 - 5) + 2(1 - v 7 ! 3 ^) 2 

PS 



In a+ + (2 - ef 



-Li 2 



aq 



+ Li 2 



-2V1 -e 2 In 



1 - a 2 + /P 2 



+ 4Vl - 5Vl - 5 - e< 



(Al) 



For the four values of A used in K° L -> £+£"7, A = 7.67, 5.67, 3.67, and 1.67 MeV the 
fractional radiative corrections are 0.0127, 0.0145, 0.0172 and 0.0223. 
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FIGURES 





(c) 

FIG. 2. Diagrams for the radiative corrections to the one-loop model are shown. A dashed line 
denotes a meson, a wavey line denotes a photon and a solid line denotes a lepton. 
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FIG. 3. Typical cut diagrams for the contributions to the absorptive part of the one-loop model 
are shown. The dot-dashed lines indicate the propagators which are put on mass shell. The sum 
of these two diagrams determines the imaginary part of the diagram in Fig. 2 (a). Figs. 2 (b) and 
2 (c) have similar cuts. 
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K -> e + e 7 
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FIG. 4. The e + e invariant mass (s) distribution for the decay K® — > e + e 7 is shown normal- 



). The solid line is the result of Ref. 



the dashed line is 



ized to the decay width To(K^ — > e + e 
the contribution of the 1/lo poles, Eq. (||), and the dot-dash line is the result for the loop model of 
Sec. III. 
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TABLES 



A (MeV) 


Modified Ref.(2) 


Loop 


Model Indep. 


7.64 


0.161 


0.154 


0.120 


5.67 


0.178 


0.172 


0.137 


3.67 


0.203 


0.197 


0.161 


1.67 


0.249 


0.244 


0.207 



TABLE I. The fractional radiative correction for iT£ —* e + e~ , as defined by Eq. (|l|), for several 
values of the cutoff A is shown. The K mass is taken to be 497.67 MeV. The second column 
is given by Eq. (^), the third column is given by the model of Sec. Ill and the fourth column by 
Eq.©. 



A (MeV) Modified Ref.(2) Loop Model Indep. 

Ym 0.0217 0.0224 0.0171 

5.67 0.0244 0.0251 0.0198 

3.67 0.0284 0.0290 0.0236 

1.67 0.0357 0.0363 0.0309 



TABLE II. Same as Table I for K\ -» n+yT 



A (MeV) Modified Ref.(2) Loop Model Indep. 

4 0.098 0.098 6.71 x 10" 2 

3 0.111 0.111 7.91 x 10~ 2 

2 0.129 0.129 9.67 x 10~ 2 

1 0.161 0.161 0.128 



TABLE III. Same as Table I for vr° e+e". We use = 135 MeV. 
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